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Abstract 

We present a study on buckling of colloidal particles, including ex- 
perimental, theoretical and numerical developments. Oil-filled thin shells 
prepared by emulsion templating show buckling in mixtures of water and 
ethanol, due to dissolution of the core in the external medium. This leads 
to conformations with a single depression, either axisymmetric or polygo- 
nal depending on the geometrical features of the shells. These conforma- 
tions could be theoretically and/or numerically reproduced in a model of 
homogeneous spherical thin shells with bending and stretching elasticity, 
submitted to an isotropic external pressure. 

Pacs : 46.70.De (Static buckling and instability), 82.70.Dd (Colloids), 89.75.Kd 
(Patterns) 

1 Introduction 

Anisotropic colloidal particles made using spheres have been the subject of var- 
ious studies in recent years. These types of colloids can be obtained in very 
different ways, for example either by deformation of spherical particles [H [2j [3], 
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or by forming clusters of them[4j El [7]. Such objects are good candidates to 
generate anisotropic colloidal crystals. Photonic bandgap (PBG) calculations 
showed that such crystals should improve the expected performances (larger 
bandgap, more convenient wavelengths) [H [8j [9]. Colloids with a hollow interior 
(spherical shells) are particularly interesting also for their mechanical properties, 
which make them potentially important for a variety of applications, such as 
drug delivery, catalysis and biotechnology, and when filled with gas as contrast 
agents for ultrasound or echographic imaging jTUl E]. Due to the relevant engi- 
neering situations as well as the biomechanics, the problem of the deformation 
of a spherical shell under external constraints has been recently investigated, 
both experimentally and numerically [T2l fT3l fT4l fT5l [TBI fT7l fT8] . 

In this study, we present anisotropic colloids obtained by buckling of spheri- 
cal shells. The buckling was induced by dissolving or evaporating the solvent 
enclosed in the slightly porous shells. This causes a stress comparable to an 
isotropic pressure on a spherical airproof shell[19j. After postbuckling, the col- 
loidal shells show bowl-like conformations, either with axisymmetric or polyg- 
onal symmetry. We observed them using both transmission optical microscopy 
and transmission electron microscopy (TEM), and compared them with config- 
urations obtained from Surface Evolver simulations, using a model of homoge- 
neous elastic (bending and in-plane stretching) spherical shells submitted to an 
isotropic external pressure. 

2 Buckling under evaporation and in solution 
2.1 Methods 

Spherical colloidal shells were prepared following previous work by Zoldesi and 
co-authors [20; 21J. They vigorously mixed dimethyldiethoxysilane (DMDES) 
and an aqueous solution of ammonia (NH3), providing a very monodisperse 
and stable oil-in-water emulsion with droplets of micrometric size. The oil 
consists of low molecular polydimethylsiloxane (PDMS) oligomers. By adding 
tetraethoxysilane (TEOS), a solid shell forms at the surface of the droplets, 
consisting mainly of PDMS with average oligomer length 4, crossed-linked with 
hydrolyzed TEOS units. As shown in ref [21], the shells are porous and then 
allow small molecules to pass through. As the low molecular PDMS can be 
dissolved in ethanol [22], adding ethanol (an equal volume in the present study) 
to such an aqueous suspension of shells filled with oil, leads to dissolution of the 
encapsulated oil into the external medium. For thick enough shells ("capsules" 
and spheres in the nomenclature of ref |21j, i.e. shell thickness over 100 nm), this 
leads to a suspension of solvent-filled spherical particles. These particles are then 
sedimented by centrifugation and redispersed in ethanol. They may be dried 
afterwards and observed by electron microscopy, exhibiting shape modifications 
("buckling under evaporation") or not OEI], depending on the ratio of their 
shell thickness to radius [T5] . 
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Thinner shells, however, will already buckle in solution when ethanol is added. 
Thin shells consist of much larger PDMS oil droplets coated with a very thin 
solid organosilica layer. They were prepared using the same procedure as de- 
scribed in ref.[20] for "microballoon" particles. In order to obtain bigger parti- 
cles, we increased the concentration of DMDES and ammonia up to 5% v/v, and 
the droplets were allowed to grow for three days before the encapsulation step. 
This resulted in somewhat more polydisperse particles with diameters between 3 
and 6 /im. Since this small polydispersity prevented from determining the shell 
thickness through static light scattering, we had to use alternative techniques 
to estimate it. A range of 5-20 nm was proposed [20] by considering that all 
the TEOS forms a dense silica shell, but careful observations of transmission 
electronic microscopy (TEM) pictures suggest 10-40 nm, which is likely since 
PDMS oligomers are known to co-polymerize with TEOS, hence contributing 
to the shell thickness. 

2.2 Experimental results. 

In order to resolve the post-buckling structure by transmission optical mi- 
croscopy, we used batches of larger particles with radii R between 2 and 3 
/im. Once buckled in solution through addition of ethanol as described in the 
previous section, these particles hold a single depression with significant vol- 
ume compared to that of the initial sphere. Furthermore, in some cases the 
depression is not axisymmetric anymore. In one batch (fig. [1]), we did ob- 
serve apparently identical objects, all with an elongated depression that gives 
the whole object a coffee-bean shape. This type of conformations was reported 
in literature for red blood cells |23j and for polystyrene shells previously filled 
with organic compounds and evaporated in air [24l [25], and was reproduced by 
numerical simulations |26j. 

In other cases, the depression presented a polygonal aspect due to regularly 
spaced radial wrinkles, in a number varying from 4 to at least 8 (fig. [2]), 7-8 
being the upper limit that we could still distinguish and count for these shell 
sizes with optical microscopy. 

The obtainment of such structures in solution from synthetic colloids is a total 
novelty. Besides, shapes with a depression presenting a 3-fold symmetry were 
previously observed in red blood cells [23} [26| [24] or in dried polymer particles 
[24} l25 j but we could not find in the literature observations concerning a higher 
number of wrinkles, on any system. To our knowledge, theoretical predictions 
leading to such wrinkles do not exist either, except when generated by a point 
[271 128] or a flat [29] load, which is not the case here. 

It is interesting to find out whether some shell heterogeneity has to be invoked 
to explain such a non-symmetry, or if a model of elastic homogeneous spherical 
shells is sufficient to recover these nontrivial shapes. For this purpose, we are 
interested in a model of thin spherical shells with bending and 2D stretch elas- 
ticity, that we derived both analytically and numerically. Such an approach was 
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Figure 1: Typical transmission optical microscopy images of particles from a 
suspension of colloidal spheres filled with oil, in a mixture of water and ethanol. 
Three different objects are displayed under three different views; they all show 
the same coffee-bean buckled shape, with an elongated depression. Scale bar 4 
jim. 




Figure 2: Spherical shells buckled in solution, showing 4 to 8 wrinkles (a-e). 
Each sub-figure shows different transmission optical microscopy views of the 
same object. Scale bar 5/im, except subflgure c: 2\im. 
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likely to give hints on features hardly accessible by experiments, such as, here, 
shell thickness, or the successive steps that lead to the final (and observable) 
shape adopted by the shells. 



3 Theoretical part: thin plate elasticity 

The elastic energy stored in the deformation of a thin sheet of an isotropic and 
homogeneous material may be split into a bending and a stretching part [30] . 
and both can be written in terms of surface elasticity: 



i 1 * •■ 1 

shell surface 



(-«(c-co) + -eijKijkie k i)dS (1) 



where k is the bending constant, cq the spontaneous curvature of the shell 
(which is zero for an unstress flat sheet, but 1/R for a spherical shell without 
constraints) |31j, and and respectively the two-dimensional strain and 

elasticity tensors. The nonzero terms of the two-dimensional elasticity tensor 
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(with A the equivalent two-dimensional Young modulus and v the equiva- 
lent two-dimensional Poisson ratio) |33j. The stretch elasticity term can thus 

be rewritten as 2 (i+iy) Tr(e 2 ) + v{ \I_l! • We ignored the gaussian curvature 
term|32j since, according to the Gauss-Bonnet theorem, its integral depends 
only on the topology for a closed surface. 

To establish the link between these two-dimensional elastic parameters and the 
three-dimensional features (sheet thickness d, Young modulus^ and Poisson ra- 
tio a of the bulk material), we follow Landau's approach[33j. For what concerns 
bending, the distribution of stress in a thin plate under flexion at equilibrium 
leads: 



(2) 



12(1 -a 2 ) 

The 2D Young's modulus A was not taken equal to Ed because this would cor- 
respond to a "planar" deformation in Landau's terminology, i.e. a deformation 
of plates at constant thickness, which seems not adequate here. The homoge- 
neous "longitudinal" deformation, i.e. without constraints in the perpendicular 
direction, is indeed more adapted to what happens in the shells. This leads: 
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It is worthwhile to notice that, as the three-dimensional Poisson ratio a has a 
maximum value of 0.5 (incompressible materials), the two-dimensional Poisson 
ratio v of such free plates cannot exceed 1/3. This point was neglected in previ- 
ous work [34], and can become of some importance if one wants to make a link 
between the parameters choosen for simulation and the geometrical properties 
of the shells. 

As the surface integral scales like i? 2 , the dimensionless Foppl-von Karman 
number 7 = is likely to drive the succession of configurations resulting 
from the balance between bending and stretching. In this model of thin shell of 
elastic isotropic material, we then expect: 



7 = 12(1- 




It is interesting to note that this model predicts conformations to be independent 
of E, and to finally depend only on the relative thickness d/R and the Poisson 
ratio. 

In such a model, we can calculate the elastic energy of an initially unstrained 
spherical surface which inner volume decreases by AV its initial value V, in two 
conformations: when the sphere remains spherical, and following reference [35] 
when an axisymmetric depression is created by inverting a spherical cap: 
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where j| can be expressed as a function of 7 and v through equation 03 Param- 
eter a is the half- angle of the revolution cone apexed at the sphere center, and 
in which the axisymmetric depression inscribes. This half- angle relates to the 
relative volume variation through: 



= i(l — cos a) 2 (2 + cos a) 



(6) 



In the limit of very thin shells and small volume variations, one can show that 



u, 



axisym 



12 l/4 
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Or, as a function of the relative volume variation: 
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In this limit, V r sphere — Uaxisym would then happen for relative volume varia- 
tions: 

— oc 7 5 (7) 

which provides a scaling law for the "sphere towards capsule" buckling. 

These theoretical calculations will be compared to simulation results in para- 
graph 14.41 

4 Simulations 
4.1 Modus operandi 

The simulated configurations presented hereafter were obtained using the free 
software Surface Evolver [36], in which the elastic energy given by equation 
(PQ) is minimized in the space of conformations. The stretch energy term is 
in fact calculated using the Cauchy-Green strain tensor, which is accurate for 
describing deformations of larger amplitude. The minimization was performed 
by alternating gradient, conjugate gradient and hessian methods. Stochasticity 
was introduced by jiggling the position of the vertices at the beginning of each 
minimization (i.e. at each volume step when the volume is decreased by steps). 
We tested that the number of vertices is high enough to avoid an influence of 
the mesh on the conformations. Furthermore, we checked that the symmetry of 
the mesh has no influence on the position of the wrinkles, by using an isotropic 
randomized mesh. Such an approach was initiated by Tsapis et al [34] . 

We explored a discrete range of 7, and restricted our simulations to the case of 
an incompressible material, i.e. v = 1/3 (a = 1/2). 

Given these elastic parameters, a first set of minimization was performed through 
stepwise decrease of the inner volume of an initially spherical surface, and stress 
free (i.e. Co = 1/R), with a minimization at each volume step. This leads to 
roughly isotropic structures with depressions regularly spaced on the surface, 
such as in reference [34] or on figure [3]-a. The number of depressions is found to 
increase with 7. 

Different shapes of much lower energy (for the same elastic parameters) could be 
obtained through more sophisticated minimizations. This was done by reversibly 
acting on the spontaneous curvature Co of the shells. Since the shells are formed 
by templating on the oil droplets, one can assume that they are unstrained in 
their initial state and cq is expected to be cq = 1/R. But when cq is changed 



Figure 3: Conformations numerically obtained with Surface Evolver for 7=743, 
z/=0.333 (which corresponds to a shell of relative thickness d/R=0.116 of an 
incompressible material with a=0.5),^=0.171, cq — l/R, 3594 vertices. Left: 
simulation performed by directly minimizing at Co = l/R for each decreasing 
volume step: N=17 depressions. Right : simulation performed first with Co = 0, 
then Co = l/R. The elastic energy is 3.10 times larger for the "potato" (left) 
than for the capsule (right). 

to zero , conformations qualitatively different could be reached: one obtains 
"capsules" with a single axisymmetric depression. Minimizing again with Co 
back to l/R preserves this capsule conformation, with an energy lower than 
the potato shape, as exemplified on figure [3]-b. Temporarily imposing a zero 
spontaneous curvature is likely to lower the energy barrier for the merging of 
two different depressions at the surface of the sphere, since merging happens 
through flattening of the high positive curvature ridge that separates the two 
depressions. This trick apparently helps to get out of some local minimas in 
which the simulated conformations are easily quenched, as is quite usual in 
buckling problems. 

In the following, such a zero curvature cycle was systematically performed at 
each volume step, in order to facilitate conformation changes. 

4.2 Results 

We performed simulations of an elastic closed surface, initially spherical and 
unstrained, which inner volume is decreased by volume steps dV = 0.0190 x Vo, 
Vq being the volume of the initial sphere. In all the experiments, Vo and n were 
kept unchanged and 7 was varied from 271 to 29160 by changing A. 

Simulations are stopped when the surface interpenetrates, which happens for 
inner volumes of the order 7-11% the initial volume Vo. 

Up to 7=583, the volume decrease causes a buckling toward the "capsule" con- 
formation, i.e. with a single axisymmetric depression, until the surface inter- 
penetrates. 
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From 7=933, axisymmetric capsules undergo a second transition when the vol- 
ume goes on decreasing, toward a non-axisymmetric conformation. The onset 
of this second transition is harder to detect since the depression, except for the 
highest values of 7, only slowly evolves toward a polygonal shape, and then 
the corners of the polygon turn into wrinkles (limits of the depression concave 
between the apices) . Quantitative data concerning buckling thresholds will nev- 
ertheless be given in the paragraph ^. 41 We could this way obtain such wrinkled 
bowls with 5 wrinkles or more. Some trials with re-increasing volume from a 
wrinkled state also provided conformations with 3 and 4 wrinkles. 

More quantitatively, figure [H provides the number W of wrinkles observed 
throughout a stepwise volume decrease for simulations with different 7 values. 
One can notice that there is some fluctuation on W. of one, or even two, units. 
In fact, performing decreasing volume simulations with the same parameters 
does not always lead to the same number of wrinkles. We could get in these 
situations an order of magnitude for the energy difference between two close 
conformation (i.e. W = ±1): it can be as small as a few tenths of percent. 
The energy of a conformation is not strictly determining its occurence: the 
path followed in the space of conformations has some importance. This is why 
we restricted most of our study to a sequence of minimization that is likely to 
reproduce the experimental situation, i.e. decreasing the volume step by step. 

Figure [4] shows a weak tendency for the number of wrinkles W to increase with 
the relative volume variation once the buckling has occured, in the same 
way that was observed in macroscopic indentation experiment s|28j. 

More obvious is the variation of W with 7. In order to precise a variation of a 
few units on a discrete quantity, we averaged W on a range of ^f- where the 
conformation holds wrinkles for all the values of 7, i.e. between 0.53 and 
0.76 (figure EJ. This puts into evidence an increase of W with increasing 7. 
Wrinkles being more numerous with decreasing j| goes in the sense of intuition: 
a thinner plate folds more easily, and hence makes more folding patterns. 

4.3 Comparison between simulations and experiments 

The inner volume of the conformations shown in figure [2] is not easy to determine 
precisely. Nevertheless, an important experimental observation is that buckling 
deformations never relaxed back toward the initial spherical shape after com- 
plete evaporation or complete dissolution of the inner oil, yet it means that the 
last water/air or oil/water interfaces, that were pulling the shell inwards, have 
disappeared. It is then likely that the shrinking brings opposite surfaces close 
enough to one another to be sensitive to Van der Waals attraction, which would 
stabilize the buckled conformation against elastic shape recovery after vanish- 
ing of the capillary forces. This hypothesis, of initially opposite parts of the 
shells that contact in the conformations experimentally obtained, seems to be 
confirmed by confocal pictures of buckled shells (figure [6]) . 
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Figure 4: Typical evolution of the number W of wrinkles held by the single 
depression after secondary buckling (simulations performed for 7=2333 (light 
grey, squares), 7465 (medium grey, diamonds), 15163 (dark grey, downwards 
pointing triangles) and 29160 (black, upwards pointing triangles; z^=0.333). 
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Figure 5: Number W of wrinkles held by the single depression after secondary 
buckling, averaged between = 0.53 and = 0.76 (error bar is the standard 
deviation), as a function of the dimensionless Foppl-von Karman number 7 = 

AR 2 
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Figure 6: Three different views of a buckled spherical shell labelled with RITC, 
in confocal fluorescence microscopy (left) and in transmission microscopy (right). 
The shell is clearly self-contacting at its convexe part. 
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Figure 7: (a) Capsule obtained from evaporation in air of a spherical shell of 
mean radius 870 nm, shell thickness 150 nm (TEM image), (b) Simulation 
obtained for 7=271 and ^=0.333 (hence equivalent to | = 0.172), ^=0.854, 
3594 vertices, (c) Simulation: 7=2333 and ^=0.333 (equivalent to j| = 0.0586), 
^=0.854, 3594 vertices. This conformation is to be compared with figure Efc 
(d) Simulation:7=20995 and ^=0.333(equivalent to | = 0.0195), ^=0.854, 
3594 vertices. This conformation is to be compared with figure^. 




Figure 8: (a) Shell enclosing oil, evaporated in air (d/R w 0.012). Transmission 
optical microscopy, size 17 fim x fim. (b) Simulation obtained for 7=2916 and 
z/=0.333 (equivalent to j| =0.064), by increasing the volume from ^=0.474 to 
^=0.209, 3659 vertices. 



In order to compare the shapes obtained through simulations (with decreasing 
volume) with the experimental ones, we thus focused on the shapes obtained 
just before self-contact. 

Figures 0-a and [7]-b shows that we could accurately reproduce the shape of 
axisymmetric capsules. In the simulation displayed here, we took 7=271 and 
^=0.333, which corresponds to d/R=0.172. 

For wrinkled bowls (examples displayed on figures 0-c and 0-d) , the conforma- 
tions obtained just before interpenetration are also very similar to shapes ob- 
served experimentally (fig. [2]-c and[2]-e). Like in the experiments, wrinkles do ap- 
pear for thinner shells. Furthermore, the parameters for which simulations pro- 



vide wrinkles are consistent with the shells' characteristics: 




to 0.08, to be compared to the experimental value 4=0.003 to 0.04. 
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Other experimental structures could also be reproduced by simulations, such as 
the one displayed on figured^. This latter was obtained through evaporation in 
air of a shell still containing some of its inner oil (note: in this case the ethanol 
was added in the aqueous solution later than usual, which possibly allowed 
polymerization of longer chains [21] that cannot be dissolved by ethanol [37]). 
Here the inner volume obviously does not correspond to shell self -contact. In 
such a process, a shrinkage of the shell itself when the water has fully evaporated 
can be invoked to explain a behaviour comparable to a volume increase, such 
as in the simulation of figure The conformation, in this case, is stabilized by 
oil- air interfaces. 

All these results show that there is no need to invoke shell heterogeneity to 
explain the shapes experimentally observed: bending and in-plane stretching 
elasticity suffices. Next section provides more quantitative insights on the sim- 
ulation of elastic buckling. 

4.4 Quantitative comparison between simulations and elas- 
tic theoretical calculations. 

The software Surface Evolver used to perform simulations provide the elastic 
energy of each conformation. For the two conformations "sphere" and "axisym- 
metric depression", we compared this elastic energy with the theoretical expres- 
sions U sp here (equation [4)) and U ax i syrn (equations [5] and [6]) . The numerical data 
are very well fitted by the theory, as shown on figure [U One sees that this first 
buckling from a spherical shape to a conformation with a single axisymmetric 
depression occurs with some hysteresis, i.e. for volume variations higher than 
the one corresponding to U sphere = U ax i sym . 

Figure fTUl presents the buckling occurences as a function of 7. The first buckling 
is determined without ambiguity, as it is obvious from figure The occurence 
of the second buckling, from an axisymmetric capsule to a "wrinkled bowl" 
conformation, is less easy to detect since it corresponds neither to a discontinuity 
nor a singularity in energy. We detected in fact two caracteristic values for , 
by visual observation of the conformations: the first one corresponds to the loss 
of the axisymmetry, when the rim of the depression becomes polygonal. Then 
the apices of the polygon becomes sharper (they tend to form the extremity of 
a c?-cone [28, 38J), and the inner part of the rim becomes convex between two 
successive apices: at this point we consider that the conformation holds wrinkles, 
and this second "threshold" is recorded. Figure [10] shows that both values are 
quite close and present the same power-law in 7 -1 . Extrapolation intercepts 
with ^ = 1 at a value 7 C =850, which is consistent with our simulations showing 
that secondary buckling appears for 7 between 583 and 933 ( j| w 0.1). This can 
be compared with calculations of ref [27] that forecast a threshold 7 of 1345 for 
the apparition of wrinkles on a clamped cap submitted to concentrated load. 

The first buckling (sphere toward axisymmetric capsule) happens for thresholds 
values of the relative volume variation that vary in a power-law with the Foppl- 
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Figure 9: Elastic energy, adimensioned by k, of a shrinking shell initially 
unstrained (co = Black squares: result of Surface Evolver simula- 

tions, with 7=2916 and ^=0.333. Continuous line: normalized elastic energy 
U S ph/K = jU sp h/(AR 2 ) of a shrinking spherical shell, calculated with same 7 
and v according to equation [U Interrupted line: adimensioned elastic energy 
Uaxisym/K of a capsule (deformation with a single axisymmetric depression), 
calculated with same7 and v according to equations [5] and [U 
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Figure 10: Primary buckling (black squares) from sphere to capsule: relative 
volume at which the axisymmetric depression appears in the simulations, for 
different 7. Interpolating straight line : = 2.8 7 -0,55 . Secondary buckling: 
polygonization of the circular rim (black upwards pointing triangles) and ap- 
parition of wrinkles (gray downwards pointing triangles) . Interpolating straight 
line: ^ = 850 7- 1 . 
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von Karman parameter: ex j~ - 55 . Despite the slight hysteresis in the 

primary buckling, this exponent is very close to the —0.6 theoretically proposed 
in equation [3 

5 Discussion. 

Experimental and numerical results showed that wrinkled bowls are preferen- 
tially observed when a very thin spherical shells lowers its volume. This con- 
formation is quite different from the structures (typically discocytes or stoma- 
tocytes) usually obtained for vesicles, where the in-plane elasticity is liquid-like 
(related parameter: 2D compressibility) |39j. Wrinkles evidences the 2D-solid 
nature of the shells, since it is needed to accomodate the surface of one hemi- 
sphere within the other hemisphere, without an excessive cost in stretch energy. 

It is interesting to note that a structure presenting 3 wrinkles had been obtained 
by Lim et al. [26] in simulations of red blood cells where the elastic properties 
of respectively the cytoskeleton and an homogeneous asymmetrical phospho- 
lipid bilayer were included in a similar numerical model with elastic bending, 
spontaneous curvature and stretching|40j. But the bending/stretching ratio in 
these biological objects, where bending and stretching have different origin, was 
higher than the range of similar values for a thin shell of isotropic material. This 
probably prevented these authors from obtaining shapes with more wrinkles. 

Besides, the simulations presented here do not necessarily provide the energies 
of lowest configuration. As an example, we could, by following another path 
in the phase diagram of elastic and geometric parameters, obtain a totally new 
conformation of much lower energy than the wrinkled bowls (figure fTTj) . But this 
conformation very likely corresponds to an energy trough too narrow to have 
been encountered in our experimental situation. Anyway, we are not looking 
for equilibrium configurations: we are trying to understand what really happens 
when a colloidal shell shrinks. It is well-known that many buckling conforma- 
tions can be quenched in non-absolute energy minima. Our study, putting into 
evidence qualitative as well as quantitative convergences between experiments, 
theory and simulations, strongly suggests that our simulations with a progressive 
decrease of the inner volume can reproduce the path followed by the buckling of 
real shells. The shapes observed are compatible both with self-contact, which 
would explain their stabilization, and with shell homogeneity. Besides, the con- 
formation (and furthermore the number of wrinkles) gives an indication on the 
shell relative thickness range. 

Structures comparable to wrinkled bowls have already been observed experimen- 
tally on millimetric half-spheres submitted to a localized[28 ; 29J or a planar|41j 
load, theoretically forecasted [271 [29], or numerically obtained by simulation of a 
sphere adhering on a flat surface [42]. but here we did put in evidence that such 
structures can also be obtained with an isotropic force distribution. 
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Figure 11: Simulation with 7=15163, ^=0.333 and ^-=0.645. Such unusual 
conformation can be obtained from a capsule at lower 7; its energy is 1.42 times 
lower than the configuration with wrinkles (W=7) obtained through progressive 
volume decrease for similar parameters. 

6 Conclusion 

Non-trivial buckled shapes were obtained by evaporating or dissolving the sol- 
vent enclosed in porous colloidal shells. 

We have shown that the deformations of such objects are consistent with a 
model of homogeneous thin spherical shells with bending and in-plane stretch- 
ing elasticity submitted to an isotropic external pressure. The numerical sim- 
ulations showed that a primary buckling leading to capsules (holding a single 
axisymmetric depression) can be followed by a secondary buckling where the 
depression wrinkles. This happens for decreasing volume variations when the 
relative thickness of the shell is reduced, and the number of wrinkles concomi- 
tantly increases. Simulations and experiments qualitatively and quantitatively 
confirm each other. 

These new results suggest that evaporation or dissolution of inner solvent is a 
promising way to obtain, from a monodisperse enough population of colloids, 
a monodisperse suspension of anisotropic objects with geometric parameters 
tunable through the characteristics of the initial spherical shell. 
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